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Abstract 



For a one-dimensional discrete Schrodinger operator with a weakly coupled potential 
given by a strongly mixing dynamical system with power law decay of correlations, we 
derive for all energies including the band edges and the band center a perturbative formula 
for the Lyapunov exponent. Under adequate hypothesis, this shows that the Lyapunov 



, exponent is positive on the whole spectrum. This in turn implies that the Hausdorff 

dimension of the spectral measure is zero and that the associated quantum dynamics 
grows at most logarithmically in time. 

o 

^. 

Ph 1 Introduction 

Let E be a topological space and Q = E^ the associated Tychonov product space. Furthermore 
let P be a propability measure on Q which is invariant and ergodic w.r.t. the left shift S : Q —>■ 
Q. Now given a measurable real-valued function V on Q and a coupling constant A > 0, one 
can associate an ergodic family of Jacobi matrices (Hx^^en (also called discrete Schrodinger 
! operators) each acting on £ 2 (Z): 

H x , w |n) = \n + 1) + A V(S n cu) \n) + \n-l), (1) 

where |n) is the Dirac notation for the state in £ 2 (Z) localized at site n G Z. If P = p® z is 
a product measure of a compactly supported probability measure p on E so that the random 
variables of the sequence (V(S n uj)) n& i of potential values are independent, the model exhibits 
so-called Anderson localization, namely the spectrum of H\ iUJ is P-almost surely pure-point 
with exponentially localized eigenstates |PFj and the induced quantum dynamics is bounded in 
time (in the precise sense given below) . The question considered in this work (and many others, 
see the reviews | Jit ^ IDamj and references therein) concerns the spectral properties as well as the 
quantum dynamics in situations where P is not a product measure so that the random variables 
(V(S n uj)) n( zi are correlated. This situation typically arises when the dynamical system (fl, S, P) 
is the symbolic dynamics associated to a (possibly weakly) hyperbolic discrete time dynamics; 
then E is the Markov partition. If now the correlations of the potential decay sufficiently 
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fast, then one expects that the model is still in the regime of Anderson localization. Here we 
complement on the prior work |CSt IBS] and prove that this holds at least in a weak sense when 
the correlations satisfy a power law decay. 

The proof of localization for these models is based on the positivity of the Lyapunov expo- 
nent. This positivity can either be established by Kotani theory [Dam], a version of Fursten- 
berg's theorem for correlated random matrices (work by Avila and Damanik cited in [Dam]) 
or by a perturbative calculation (for small A) of the Lyapunov exponent. This latter calcula- 
tion was first done by Chulaevski and Spencer [CSj by carrying over the argument of Thouless 
|Tho] . in a version given by Pastur and Figotin |PFj . to the case of correlated potential values. 
The resulting formula is recalled in Section [2j Based on this result, Bourgain and Schlag then 
proved localization |BSj . The only flaw left is that in |CSj (and actually already in |PFj ) not 
all energies could be dealt with, but the band center and the band edges were spared out. Here 
we show how the techniques of our prior works on anomalies and band edges |Sch[ ISSj combine 
with those of |CS] to rigorously control the perturbation theory for the Lyapunov exponent also 
at these energies. Instead of repeating the rather complicated proofs of |BSj . we then adapt 
to the case of correlated potentials the elementary and short argument of |JSj showing that 
positivity of the Lyapunov exponents implies at most logarithmic growth of quantum dynamics 
and hence, by Guarneri's inequality |Guaj . zero Hausdorff dimension of the spectral measures. 
Even though this is a weaker localization result than pure-point spectrum with exponential lo- 
calized eigenf unctions, it proves the behavior which is stable under perturbation and we hence 
consider, as argued in [JSJ, that it already captures the physically relevant effect. In the next 
section, the results and the precise hypothesis are described and discussed in detail. The other 
sections contain the proofs. 

Acknowledgment: This work was supported by the DFG. 

2 Set-up and results 

In order to fix terms and notations, we have to begin by reviewing some basic definitions of 
symbolic dynamics and strong mixing [Bowl IPPj . Let S be a countable set furnished with 
the discrete topology. We designate a reference element G X. For any subset / C Z and 
u = (cn)nez £ let us define 

iti(uj) = (a n )n<zz , °n = for n g I , a n = a n for n G I . 

For a bounded, measurable function g : Q — > V into a real, normed vector space (V, || . ||), the 
variation on I is defined by 

Var 7 (#) = sup \\g(u) - g(u')\\ ■ 

717 M=17 (^') 

Then g is called quasi-local with rate < r < 1 if and only if there exists a constant C = C(g) 
such that, for any m,n > 1, 

Var hm>n] (<7) < C(g)r mAn } m An = min{m,n} (2) 
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The set of all quasi- local functions with rate r is denoted by Q r (Y). 



Next let us state precisely the strong mixing hypothesis used in this work. For m < n and 



ak G X with m < k < n, the associated cylinder set is denoted by A 



m,n 



A-r. 



{u = (<Jk)kei \ <T k = a ki m < k < n). Then the invariant measure P on the shift space (ft, Z) 
is said to satisfy a power law -^-mixing |Braj with exponent a > if there is a constant C > 



such that for all k < I < m < n and all A^i, A m ^ n , one has 



V{A Kl n An,n) - P(A k)l )P(A 



m,n ) 



< CP(A k>l )P(A 



m,nj 



\m 



l\ 



Equivalently, for any vr^/j-measurable function g\ and 7T[ mjn ]-measurable function g% holds 
E(g l92 ) - E( gi )E(g 2 ) < C E(\ 9l \) E(\g 2 \) \m - l\~ a , 



(3) 



(4) 



where k < I < m < n and C as above. This also implies ergodicity. Examples when (J3J) holds are 
given in Remark 1 and 2 below, after the main results are stated. Averages over uo w.r.t. P are 
denoted by E, or also by E^ if the dependence on uj is retained in the integrand. Furthermore, 
the set of centered quasi-local functions will be denoted by Q°(V) = {g G Q r (V) | E(g) = 0}. 

Throughout we suppose that the potential in (JTJ is given by a centered real-valued quasi- 
local function V G Q°(R). It is well-known and verified in Lemma H] that (j4j) implies the decay 
of correlations \E(V(u))V(S n u)))\ < C\n\~ a for some constant C. For a > 1, one can hence 
define its (positive) spectral density Dyik) at k G [0,27r): 

N-l 



Dv(k) = ^e !fa E4y( W )y(S" W )) 



lim 4; E u 

N ^oo N 



\J2 e lkn V{S n u) 



n=0 



As final preparation let us recall the definition of the Lyapunov exponent 'jx(E) at energy 
E G C associated to (JT|). If the transfer matrices are defined by 

'E - XV {uo) -V 
1 



T E 

2 X,uj 



Q r (SL(2,R)) 



(5) 



then 



A' 



lim — E u 

TV^oo N 



n=l 



E 

X,S n uj 



The main result of Chulaevski and Spencer |CS] is that for a > 2 and at an energy E 
2 cos(A;) in the spectrum [—2, 2] of the discrete Laplacian, but away from the band edges E 
—2, 2 and the band center E = such that the distance d{k) of k from mod | is positive: 



lx{E) 



A 



D v (k) 



+ O 



3a + 2 

A a + 2 



(6) 



8 sin 2 (fc) ' w V d(A;) 

As we need to build up the whole formalism anyway, the main element of the proof of (EI) is 
reproduced in Section El As indicated, the control of the error terms breaks down at the band 
edges and the band center. Our first result provides perturbative formulas for the Lyapunov 
exponent at these energies, generalizing respectively our prior results for independent potential 
values |Scht 
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Theorem 1 Assume a > 2, Dy(0) > and Dv(tt) > (the latter is only needed for (i)). 

(i) The Lyapunov exponent near the band center E = is given by 

7A (eA 2 ) = A 2 ^^ f d0p e (0)(l+cos(40)) + O(A^), (7) 
° Jo 

where p € is a n-periodic smooth probability density. 

(ii) Up to errors of order 0(\ 3a + 6 ), the Lyapunov exponent near the upper band edge E = 2 is 
given by 

lx (2 + e\l) = A§ j &9p e (9) sin(2#) + ^^ J d9p t (9) (l + 2 cos(2#) + cos(4#)) 

(8) 

where p e is a n -periodic smooth probability density written out explicitly in Section [HI The same 
formula holds at the lower band edge E = —2. 

The formulas (jSJ), (ED) and (jSJ combined allow to study the Lyapunov exponents at all 
energies [—2, 2]. In order to assure positivity for A > 0, one first has to check that the spectral 
density is positive (cf. Remark [3] below) and then prove that the integrals appearing in ((Tj) and 
(JHJ) are positive. This is immediate for (|7j). For (JHJ we could not produce an analytic proof, 
but, given the explicit formula (1421) for p e , one can readily do a numerical evaluation. 

Nevertheless, the three formulas are not yet sufficient to prove uniform positivity of the 
Lyapunov exponent on the whole spectrum for some fixed small, but positive value of A. Indeed, 
for once the non-random spectrum a(H\^) may (and typically will) fill the whole interval 
[—2 — A|| V||oo, 2 + A||V||oo] where ||V||oo = P — esssup |V(u;)| (use approximate eigenfunctions 
as Weyl sequences in order to show this). For an energy 2 + eA, e > 0, the asymptotics (JSJ) then 
says nothing. However, one can combine the techniques of this work with those of [SS] in order 
to prove, as in the case of independent potential values [SS[ Section 8], 



7 A (2 + eA") = Ve~X^ + 0(X 1 ~ R ^ A^" 1 , A^ w) , e >0, (9) 



where | < r] < | is such that the error terms are of lower order than A 2 (in particular, 77 = 1 
is allowed for a sufficiently large). Moreover, the formulas ([6]) and (jSJ) do not imply positivity 
of the Lyapunov exponent at a fixed A for all energies in [2 — A, 2) because the error term in 
(jSJ) explodes as one approaches the band edge. However, once again one can transpose [SS, 
Section 8] to the case of a strongly mixing potential: 

7A (2 - eA") = A 2 -" ^ffl + 0(A 4 "f , \%, , e > , (10) 

8 e 

where again | < r] < | has to assure that the error terms are subdominant. A careful anal- 
ysis now allows to show (modulo the issues discussed above) that for A sufficiently small the 
Lyapunov exponent is positive on [2 — c, 2 + c] for c > 0. We do not provide the detailed 
argument here, but do claim to have presented all the essential ingredients in order to complete 
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it. Similarly, by analyzing the Lyapunov exponent j\(e\ v ), 1 < 77 < 2, using the techniques of 
[Sch, Section 5.1] or [SSj Section 5], one can show that the Lyapunov exponent is positive near 
the band center for A sufficiently small. 

Let us now assume that uniform positivity of the Lyapunov exponent has been verified for 
all energies in the spectrum, either by the above or some other argument, and then deduce 
localization estimates from this. One standard way to quantify the spreading (derealization) 
of an initially localized wave packet |0) under the quantum mechanical time evolution e~ ttHx ^ 
is to consider the growth of (time and disorder averaged) moments of the position operator X 
on£ 2 (Z): 

roc j± 

M| = / — e-T E w (Ole^^lX^e-^- 4 10) , g>0. (11) 
Jo * 

Boundedness of M|, uniformly in time is called dynamical localization. Logarithmic growth in 
time as obtained in the following theorem is quite close to that. 

Theorem 2 Consider an ergodic family of Jacobi matrices (H\^) w( zq of the form (CQ) with a 
quasi-local potential V and an invariant measure P satisfying (jlj) with a > 0. Suppose that the 
spectrum is included in an open interval (E , E{) on which the Lyapunov exponent is uniformly 
positive: 

lx{E) > 70 > 0, Ee(E ,E 1 ). (12) 
Then for any (3 > 2 there exists a constant C ((3, q) such that 

M| < (logT)*' 3 + C(P,q) . (13) 

Furthermore, the H aus dor ff dimension of the spectral measure of H\ )U} vanishes P '-almost surely. 

The elementary proof (fitting on 4-5 pages) of ( Tl3|) is almost completely contained in [JSJ. 
It is therefore not reproduced here, but we discuss in detail in Section [9] the only step that 
has to be modified. As already indicated in the introduction, the last statement then follows 
directly from Guarneri's inequality [Gua] . 

Now follow remarks on when the hypothesis of the above theorems are satisfied. 

Remark 1 The strong mixing condition (jlj) clearly holds if P is the product measure of some 
probability measure on S, because the functions gi and #2 are then independent. The mixing 
condition also holds if P stems from a Markov process given by a stochastic kernel having 
only one invariant measure on a countable set S. Then the decay on the r.h.s. of (jlj) is 
actually exponential, with rate given by the Perron-Frobenius gap of the stochastic kernel. 
Yet more general, let us consider a hyperbolic dynamical system (X, T) (Axiom A) given by 
a map T : X — > X. Then one has a finite Markov partition S, with associated symbolic 
dynamics (Q, S) [Bow], and there is a wealth of so-called Gibbs measures associated to Holder 
continuous (i.e. quasi-local) functions which all satisfy (jlj) with an exponential mixing rate 
[Bow, Proposition 2.4]. Two standard examples of this type already cited in |CS] are the 
period doubling map and the Arnold cat maps. Moreover, if the phase space A is a manifold, 
then any differentiable real function on this manifold gives rise to a quasi-local potential under 
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the coding map. For all these examples with exponential ^-mixing, the error bounds in (J6]), 
(J7J) and ([8]) are given by the error bounds of the independent case \PF\ \Sch\ ISS] multiplied by 
log 2 (A). The error bounds in the independent case are recovered by sending a — > oo in ([6]), ([7]), 
©, © and (HUD 

Remark 2 Concrete examples of dynamical systems (X, T) having not an exponential, but 
only a power law decay in (Tj0) have only be analyzed more recently. Necessarily T is then not 
uniformly hyperbolic, but it is supposed to have only a few parabolic points. Such examples can 
be constructed even if X is an interval, but the invariant measure then has a non-normalizable 
density w.r.t. the Lebesgue measure. It is, however, possible to construct a symbolic dynamics 
over a countable alphabet £ which then has an shift-invariant propability measure P satisfying 
the strong mixing estimates ([3]) and (TJJ. Instead of producing a long citation list, we refer to the 
references in |Gou] which contains a proof of (jlj) for several concrete examples. It is precisely 
in order to deal with these cases at the verge that we bothered to work with (J4]) instead of 
exponential mixing. 

Remark 3 The positivity of the spectral density Dy(k) can for some examples be checked by 
an explicit calculation, but there are also further techniques available in order to verify this 
[Bra]. The case of Dy(0) is particularly well studied because of its importance for central limit 
theorems [PPj IGou] . For the Gibbs measures of Remark 1 and the examples of Remark 2, 
Dy(0) = holds if and only HV = voS — t> is a cocycle given by another quasi-local function 
v. By suspension, one can deal similarly with k = ir and actually any rational Jf-. 

Remark 4 The above results transpose if Z is replaced by N, namely for Q = S N furnished 
with the left shift and Hx ^ acts on £ 2 (N). As the inverse S~ l of the left shift operator is not 
defined in that case, one needs to replace in all proofs functions like goS~ n for n > by (U*) n g, 
where U* is the L 2 (£ N , P)-adjoint operator of U : g i— > g o S. 



3 Anomalies at band center and band edge 

Let us begin by recalling that the transfer matrix T X E W G SL(2,R) given in ([5]) is elliptic for an 
energy E = 2 cos(A;) G (—2, 2) and A = 0, and it can hence, to zeroth order in A, be transformed 
into a rotation. More explicitly, 

where 

sin(k) cos(k) J ' J " ^/sin(fc) V - cos ( fc ) K 

In the next section we will consider the action of the matrix (fl~4l) on the real projective line, 
which is identified with a circle. To lowest order A , this action induced by (fl4"l) is then a 
rotation on the circle. For irrational J^, there is a unique invariant measure given by the 



sin(fc) V 1 



/ cos(/c) — sin(A;)\ 1 / sin(A;) 
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Lebesgue measure. For rational = - at least Birkhoff sums of harmonics of order lower than 

° 2-7T q 

q vanish. 

At the band center k = |, the square of the transfer matrix (fl4j) (note that M = 1 here) is 
the unit matrix and one can only control the lowest order harmonic, which turns out not to be 
sufficient for the calculation of the Lyapunov exponent. It is then more convenient to consider 
directly the square of the transfer matrix 




and to group the coordinates of ui in pairs and consider Vt = Tr where S = S x S, and furnish 
it with a probability P naturally induced by P. Again the suspension (Q, P) is a shift space 
with power law mixing. However, the matrix (TT5]) is now in the form of an anomaly as discussed 
at the end of this section. 

At a band edge, e.g. E = —2 and k = ir, the basis change in (Tbil) becomes singular and one 
has a Krein collision. Nevertheless, the transfer matrix at A = can be transformed into a non 
diagonalizable Jordan normal form: 

NT ^ N ~ ~\ AV(u,)-eAl J ' l-l l) ■ 

Let us further conjugate this matrix by N\ = in order to get again an anomaly (c/. 

[SSj for a motivation of this conjugation): 

NxNT -^ AT-AT-i = - exp (a* °) + At (_° e J) + O(A)) . (16) 

Resuming, after adequate basis change and possibly regrouping of terms, one has to study 
in each of the three situations (|T^|) . (|T5l) and (|T6l) families of random matrices (T XtUJ ) \>o^en £ 
Q r (SL(2, R)) of the following form: 

= ±R k exp(A"P 1)aJ + A 2 "P 2 , W + G(A 3 ") ) (17) 

where 77 > 0, P iiW G Q r (sl(2,R)) for j = 1,2, E(P 1>W ) = and the error term C(A 3r? ) is 
uniformly bounded (i.e. the bound is ^-independent). If k = 0,7r, namely at a band center 
(|T5l) and a band edge (1T61) , such a family is said to have an anomaly of second order [Schl, ISS] . 
In the following sections, we treat general families of the form (fTTl) . and then go back to the 
explicit cases in Section [8] in order to complete the proof of Theorem [lj 
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4 Phase shift dynamics 

The bijective action St of a matrix T £ SL(2, R) on = R/7rZ = [0, n) is given by 

with an adequate choice of the sign. This defines a group action, namely Stt' — StSt>- In 
order to shorten notations, we write S\ u = St x where is of the form (TlTI) . One thus has 
S Xi U0)=9 + k + O(\). 

Given an initial angle 8q and iterating this dynamics by the left shift on Q defines a stochastic 
process 9 n (u), also simply denoted by 6 n below: 

O M = , 9 n+1 (ou) = Sx,sn u (9 n (u)) . (19) 

In order to analyze the dynamics in more detail, let us introduce for j = 1,2 the trigono- 
metric polynomials 

One then has (c/. [SSJ for details) 

2 i 

SxAO) = + k + J^^PjAO) + ^ 2n Pi^d 0Pl ^{e) + 0{\^) . (21) 
i=i 

Due to Lemma 3 of |JSS] and a telescoping argument, the Lyapunov exponent 7(A) char- 
acterizing the exponential growth of the products of matrices in the ergodic family (T\ ) gn aJ ) n > 
is given by 

^ = i im ^ E ^ E - E teg(\\Tx,s^ee ni u,)\\) , (22) 

TV — >oo iv — 
n=0 

where Ee denotes an average over the initial condition 80 w.r.t. an arbitrary continuous 
probability measure on S^. As our interest is perturbation theory for 7(A) w.r.t. A, we shall 
need the following expansions for the summands of ( |22l) (e.g. |SS] contains the algebraic proof): 



Lemma 1 Set 



a 



3,v 



(v\P j;Ui \v) , = (v\ \v) , 7^ = (v\\P j:UJ \ 2 \v) 



Then pj jL0 (9) = ^sm(aj !U} — f3j )UJ e 2ie ) . Furthermore, 

log(||T Ai ^||) = 3?e (J2X V ^e™ + \ (|^| 2 + !^e™ - fi^e**)) + O(X^) . (23) 

3=1 
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Formula (1221) and also its perturbative evaluation based on (|23|) hence leads us to consider 
sums of the type 



I N (G) = t; E « V) SOS^AM) , '(0) = Km , (24) 

iV ' iV— >oo 

n=0 

for functions £ on x of the type C/(u;, 9) = gj(u)f(9). More explicitly, the above lemma 
shows that one only needs functions of the form g(uj)e and g(u)e 4ie with g G Q r (C). For a 
7r-periodic function / G C(S'i), we also introduce 

AT-l 

= T7E^/(^ n ), /(/) = Jim /*(/), 

/V ' jV— >oo 

n=0 

This is a Birkhoff sum of the process n = # n (a>). In the sum (1211) there is, moreover, an explicit 
dependence of Q on hence let us use the term Birkhoff-like sums for the sums 1(G). 

5 From Birkhoff-like sums to Birkhoff sums 

The aim of this section is, as indicated in the title, to reduce the perturbative evaluation of the 
Birkhoff-like sums (|24|) to the evaluation of Birkhoff sums by invoking the correlation decay 

®- 

Proposition 1 Suppose a > 2. Let g G Q r (C) and f G C 2 (Sl). Define G(u,9) = g(u)f(9). 
Then 

1(G) = E(g) 1(f) + O(X^) . (25) 
IfE(g) = 0, one has the following convergent expression for the next higher order contribution: 

oo 

1(G) = A" £/(/.,) + 0{\"&) , f 3 (9) = E u (g(S^) Pl ^(9))f'(9) • (26) 

3=1 

The first lemma needed for the proof is mainly contained in |CSj . We provide a few more 
details of the proof and use the notations of this work. 

Lemma 2 One has for m, n > 1 

Var[_ m _ nin+m] (0„) < 0(r m \ r <) . 

Proof. Using equation (ITTJ|) . 

\O n+1 (uj) - 9 n+1 (u')\ < |<SA,s»^n(w)) - «SA,s« w (0 n (a? / ))| + \S\,s n w(9 n {v')) - S\ )tl3 '(0 n (u'))\ , 
one deduces 

\e n+1 (u) - e n+1 (u')\ < fsup \s' Xu {6)\) \e n (u)-e n (u')\ + ^ v \s XiS ^(e) - s XiS ^(e)\ . (27) 

\ u>,9 / e 
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Using the estimate 



x x 



x — x 



+ x' 



\x\ 



\x \\ — \\x\ 



< 



\x\ 



\x — X 



and the definition of S XuJ , it follows 



This implies 



< 



IT; 



|T\,o; - T\,o/|| < 2 ^SUp ||TjJ,||^ \\T XjU) - T XyU 



(2i 



where C\ is a constant and T X)U = 1 + X V Q X ^ for a matrix- valued function Q x <u that is analytic 
in X v and uniformly quasi-local for small A (i.e. the constant and rate is A-independent). 
Furthermore, one has 

sup 15^,(0)1 < 1 + C 2 A" 
for A sufficiently small and some constant C 2 . Applying this and (I28I) to (127]) one gets 

Va r/ (0 n+1 ) < (1 + C 2 A") Var 7 (0 n ) + Ci A" Var^Qv^) . 
Iterating this estimate and using Var/(0 o ) = 0, it follows that 

n 

—m—n,n+m] 

3=1 

oo 

< CiA"C3r m+1 ^(1 + C 2 A")V = 0(A"r m ) 
i=o 

for A sufficiently small. □ 
In order to state the next two lemmata, we introduce the following notation extending 

1 N ^ 

W) = ^E^£(s m +^AH), i m (G) = lim rsig). 

'« ' TV— >oo 



n=0 



Lemma 3 Let <7i, #2 £ Qr(C) and / G C^S^). Furthermore let k > I > and m > 1. T/ien 

E4^(5 3m+fe+ "^)^(5 3m+i+ ^)/(0„(^))) = E(/(0„,))E(( ?1 oS fc -^ 2 ) + G(m- a ), (29) 
uniformly in k, I and n. This implies, for Q(oJ, 9) = gi(S k uj)g2(S l uo)f{9), 

I 3m (Q) = E(( 9l oS k - l )g 2 )l(f) + 0(m" a ) . (30) 
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Proof. By Lemma [2] and because / is Lipshitz-continuous, one has uniformly in n 

\f(e n (ou))-f(e n (7r hm - n , n+m] (Lu)))\ < 0(A"r m ). 

As g\ and g 2 are quasi-local and therefore bounded, one also deduces uniformly in k, n and I 

9l (S k+n+3m uj) g 2 (S l+n+3m uj) - {( 9l oS k )(g 2 oS l ))oS n+3m on [n+2mtn+k+Am] (uj)\ < 0(r m ) . 

Let us denote the two functions inside the modulus by g and g respectively. Similarly denote 
/ o 9 n o 7T[_ n _ mn+m ] by /. Now consider E(g (/ o 6 n )). As the functions / and g are bounded, 
it follows from the estimates above and (J4]) that with errors of order O = 0(m~ a ) > 0(r m ) > 
0(X v r m ) (for big m and small A) in each step we get 

E(gf(9 n )) = E(gf(9 n ))+0 = E(gf) + O = E(g)E(f) + O = E(g)E(f(9 n )) + 0. 

This finishes the proof. □ 

Replacing g2(S 3m+n+l uj) by g2{S l+n u) for < I < k, one can modify the argument by 
grouping g 2 and / together. This gives the following 

Lemma 4 Let g±,g 2 £ Q r (C) and let f 6 C7 1 (S' 7 J:). Then one has for < / < k and m > 1 
E4^(S 3m+ " + M^(^ + ^)/(^M) = E(gi) E (J (g 2 (S l+n (u})) f(9 n (u)) + 0(m~ a ) , 
uniformly in I, k and n. This implies for Q(u, 6) = gi{S 3m+k oj)g 2 {S l u)f{6) 

1(G) = E( gi )l(g 2 (S l Lu)f(6)) + 0(m- a ) , (31) 
and leads, for f = 1 and I = 0, to 

E( 9l (S 3m+k (uj))g 2 (uj)) = E{ 9l )E{g 2 ) + 0(m" a ) . (32) 

Proof of Proposition [H By Taylor expansions and Pi,s n +iu(8n+j) — Pi,s n +^(^n) + 0(j\ v ), one 
finds 

6m— 1 

f(9 n+6m ) = f(9 n ) + X V Y1 Vi^ + aOn)f\e n ) + G(m 2 A 2 ") . 
Therefore multiplying with g o s 6m+n and averaging over u and n gives 

6m— 1 

1(G) = f 6m {G) + A" ]T f{Gj) + 0{m 2 \^) , 

3=0 

where Gj(u3, 0) = g(S &m ui)pi^ S j ui{0) f '{9) . As is a trigonometric polynomial in 9, Lemma[3] 

can be applied to each summand in order to obtain 

i em (G) = E{g)I{f) + 0(m- a ) . 
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Because the functions Qj are uniformly bounded, one has A'Eh ^i) = 0(mA"). Using 
m = A _T?T +" now proves the first part. 

Now let E(g) = 0. Again because px jW is a trigonometric polynomial, Lemma [3] gives, for 
j > 3m and fj as defined in (T2"6T) . 

Using Lemma HI one obtains for j < 3m 

i{G 3 ) = E(g)P(p li0 ,f(6)) + 0(m~ a ) = 0(m- a ). 

All together, one has 

6m— 1 

1(G) = A" J2 J (/6m-,) + 0(m 2 A 2 ',A% 1 - tt )m - a ) . 

j=3m 

Because (132]) gives 

oo oo oo 

j'=3m+l j'=3m+l j'=3m+l 

one therefore deduces 

oo 

1(g) = XnJ^Iifj) + 0(m 2 X 2r ',X r >m 1 ~ a ,wr oc ). 
Finally choosing m = A q + 2 concludes the proof. □ 



6 Oscillatory sums away from band center and edges 

As already explained in Section HI for the calculation of the Lyapunov exponent one needs to 
evaluate the Birkhoff-like sums of functions of the type Q(u,9) = g(uj)e 2t ^ 6 , j = 1,2. This 
is done in Proposition [2] below for energies away from the band center and band edge. By 
applying it to the terms appearing when (|2"3"|) is replaced in (T2"2"l) . this result allows to complete 
the proof of formula As the straight-forward algebraic calculations are carried out in detail 
e.g. in |CS| IJSS] and we present a similar calculation for the band edge in Section [8], we skip 
the details. 

Proposition 2 Let a > 2. Suppose that the lowest order rotation phase k in the dynamics 
( |2T!) satisfies d(k) = dist(fc mod|,0) > 0. Consider Qj(u, 9) = g(uj)e 2ljd with j = 1,2 and 
g e Qr(R). Then 
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// ; moreover, E(g) = 0, 




Proof. [PFIICSIUSS] The dynamics and the definition of the Birkhoff sums implies /^(e 2 *- 761 ) = 
e 2ljk I N (e 2lje ) + 0(N~\ A 77 ). This implies I(e 2lje ) = 0(d(k)~ 1 X' n ). The bound ([25} of Proposi- 
tion [1] thus implies the first statement. The formula (1331) now follows after a short calculation 
from (1261) . the identity P\^(9) = 3m(« lu — Pi )U) e 2td ) and the first statement. □ 

7 Fokker-Planck operator for drift-diffusion 

We now focus on energies for which the rotation angle k in (12T!) satisfies k mod | = so that 
the argument of Proposition [2] does not apply in order to calculate the Birkhoff sum I(e 2l °). 
For this purpose, let us introduce the bilinear form 

oo 

(gi,g 2 )a = E ul (g 1 (u;)g 2 (cu)) + 2 ^ E w (g 1 (uj)g 2 (S n uj)) , 9l ,g 2 G Q° r (R) , 

m=l 

which by f[3"2"j) is well-defined. Note that D v (0) = ( V, V)q. Let us use the notation Pj(u, 9) = 
Pj,u)(9) an d Pj — dePj- Then expressions like (pi,p[)n are functions of # on S^. 

Proposition 3 Let the family T\ jU1 be as in f[T7|) . and F G C 3 (Sl). For f G C 1 (Sl) given by 

f = ( Pl ,Pi)nF" + ((p^tih + 2B{p 2>ul ))F' , (34) 

one then has for a > 2 

i(f) = G(A"Si) . 

Proof. By a Taylor expansion, one has with errors of order O = (D(X 3ri ) 

F(S X ,M) = F(9) + j^X k ^(9)F>(9) + X 2 ^- [F\9)pU0)p[J0)+pU0)F"(9)] +0. 

k=i 

We now use this for 9 = 9 n and average over n. Because is centered and a polynomial, 
one can apply equation fl26|) of Proposition [1] to the term with power X v and fl25|) to the other 
terms. This gives 

1(F) = 1(F) + ^(/((^K)^') + l(( P uPihF") + 2/(E,( P2 , w )F'))+0 

with errors of order O = 0(X TJ a + 2 ). As the functional / is linear, resolving this equation for 
1(f) gives the desired result. □ 
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This proposition shows, that we can control error terms on Brikhoff sums for a function /, 
if / is in the image of the operator £ on functions on Si given by 

£ = (pde + q)d e , P = (PhPi)n , q = (pi,p[)n + 2E(p 2 ,u,) • (35) 

As one needs to calculate Birkhoff sums /(/) pertubatively, we are looking for some class of 
functions where lim^o 1(f) exists. For / in the image under C of C 3 (Sl), this limit is 0. Thus, 
if this map is given by the scalar product with some L 2 -function p, one has p G R,an(£) ± = 
Ker(£*), where the formal adjoint is given by 

C* = d e (d e p - q) . 

C* is a forward Kolmogorov or Fokker-Planck operator describing the drift- diffusion dynamics 
of the process 9 n on S^, and C is the associated backward Kolmogorov operator [Risj . It will be 
shown that in the situations considered here, Ker(£*) is spanned by a smooth, /^-normalized 
function p. Furthermore, the following theorem shows that / G Ker(£*) ± D C 2 (S\) turns 
out to be sufficient for finding a solution F G C 3 (S%) of the differential equation (134")) so that 
Proposition [3] actually applies. Even though contained in |SSj . let us give the proof for sake of 
completeness. 

Theorem 3 Suppose thatp(9) = for at most one angle 9 G S^. Furthermore suppose q(9) ^ 
in that case. Then the Fokker-Planck operator C* has a unique groundstate p G C°°(Sl), which 
is non-negative and normalized. Furthermore, for f G C 2 (Sl), one has 

1(f) = [ n d9p(9)f(9) + 0(A"fit) . 

Proof. Integrating the equation C*p = once gives 

(pd e + (d e p) - q)p = C , (36) 

where C is some real constant. As /(/ + c) = c + /(/) for c = (p,f), we may assume 
f d9 p(9) f(9) = once we found the normalized solution of f)36p . Proposition [3] then gives 
the bound on /(/) if one finds a solution G G C 2 (S' 7 J:) of 

(pde + q)G = f, I* &9G(9) = 0. (37) 

Jo 

First let us consider the case p > 0. Then there is no singularity and C* is elliptic. The 
groundstate p and the function G can be calculated. For some 9 set 

W{6) = Je^W)' W{9) = Jo ^ W) W(6) = Je d ^ e • (38) 

Then 

w 

p = C x — (C 2 W + 1) , G = e- w (W + C 3 ) , (39) 
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where C2 is fixed by the condition that p is 7r-periodic and G\ > is a normalization constant. 
This fixes C = C\ C2 in (|36j) . G is a solution of the first equation of (f3Tj) and for C 7^ the 
constant C3 is fixed by the condition that G is 7r-periodic. Furthermore one has 







j pf = J p{pd e + q)G = - J G (dop - q)p = -C j G(9) . (40) 



Thus G is a solution of (J57j) . If C = 4^ C 2 = 0, then w is ^-periodic as well as W which 
follows from j pf = 0. Therefore G is ^-periodic and C3 is chosen such that the integral in 
(1571) vanishes. 

Now let p{9) = for exactly one 9 G and for sake of concreteness let q{9) > which 
implies q{9) > 0. Then choose 9 G (0, 6* + 7r) in the first equation of ( 1551) . 9 = 9 in the second 
one and = 9 + 7r in the third one. As lim^g e w ^ = and lim ^ +7r e^ 9 ) = 00 in this case, 

u;, W and W are well-defined for 9 G (0, + 7r). Using de l'Hospital's rule, one can prove by 
induction (see [SS] for details) that 

f w 

p = C — W, G = e~ w W , 
V 

can both be continued to a smooth (even at 9) and 7r-periodic function. C > is again a 
normalization constant and hence equation (j4l3l) shows that G solves (137|) . □ 

Before applying this result in order to prove Theorem [TJ let us present another derivation 
of the equation C*p = 0, albeit a formal one, which shows that p is the lowest order ap- 
proximation for the assymptotic invariant measure of the process 9 n . Expanding the function 
= S Xt s N - 1 u) • • • <S\,Slu <S\,u> shows that the coefficients of 

S^(9) = 9 + yp N w {9) + -A 2 ^) + O(X^) , 

are 

N-l N-l / n-l \ N-l 

= ^Pi,s^ , q" = Yl [Pi^u + ^2pi,&u I Pl,S*u> + 2 ^2P2,S^ ■ 

n=0 n=0 \ j=0 J n=0 

An invariant measure u\,n for N steps of the dynamics 9 n on satisfies 

f v x>N {d9)f{9) = E r^xA^fiKs^M) > / e C{Sl) . (41) 
Jo Jo 

Supposing v\^(d9) = p\,N (9) d9 = p N {9) d9 + o(A°), gl]) leads to 

C nPn = 0, C% = d e (d e E({p^) 2 ) - E(^)) . 
Using the stationarity of P and the definitions of p^ and q^ , one deduces 

lim i-E((pf w ) 2 ) = p, lim ^E(g^) = q, 

where the convergences are uniform in 6 1 . This shows that ^jC^t C* weakly for A — > 00. 
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8 Application to the band center and band edge 



This section contains the proof of Theorem [H Let us first consider item (i), that is the band 
center. As described in Section [3] we have to work with the probability space = (Ex S) z 
which is isomorphic to Q by the pairing isomorphism V. Using this isomorphism and the 
potential V, which is defined on Q, let us define the two random variables on Q 

v* = V{V-\u)) = V{u) , = V{SV~\u)) = V(Su) . 

Then according to equation ffl5|) the family of matrices we have to consider is given by 



Txfi = - exp 



A ( ° + * (~ U f » 26 ) + 



/ 2 V — 2e ucjV t 



In this situation one has = i{vu> + u^)/2, (3\ t Q = — v^)/2, a 2 ^ = —it and /3 2j & = 
— Using Lemma [1] and (v — u,v — u)q — 2 Dy{n) and (v + u, v + u) q — 2 £)y(0), one 

obtains that the polynomials fl35l) are explicitly given by 

p(6) = l -D v {$) + ~D v (ir) cos 2 (2fl) , q(6) = -^D v (ir) sin(4#) - e . 

By assumption on V, one has p > uniformly on S^.. By Theorem [3] there is thus a smooth, 
positive and /^-normalized groundstate p e for the operator C* (which can readily be written 
out). Furthermore, one checks ji t & = (w? — u%)/2. Then equation (1231) . Theorem [3] and 
Proposition [T] combined with some algebra leads to (JTj) for 7 A (eA 2 ) = 57(A)- 

Now let us prove Theorem C](ii). Hence let T AjLJ = N x NT~ 2+eX N^N' 1 be the anomaly 
given in (USD- As a 1>CJ = iV{uj)/2 , fa*, = -iV(u)/2, a 2 ^ = -i(e + l)/2, and (3 2 ,u> = i(e- 1)/2, 
one deduces, using (V, V)n = Dy(0), 

p(0) = D v (0) cos 4 (#) , q(9) = -e - 1 + (1 - e) cos(2#) - 2D V (0) cos 3 (fl) sin(0) . 

By assumption on V one has p{6) > for 6 ^ |, and as g(|) = —2^0, there is a unique 
groundstate p t G C°°(S' 7 J:) by Theorem [3J Explicitly, one obtains 

Pe(6) = C j J£ exp ^___ (t ail 3 (0 - tan 3 (£) + 3etan(£) - 3etan(*))) , (42) 

where C is some normalization constant. Furthermore, one checks 7i ja; = V(uj) 2 /2 and hence 
(1231 . Proposition [1] and Theorem [3] imply (El). 



9 Bound on the quantum dynamics 

As already said above, the proof of Theorem [2] follows exactly the proof of Theorem 1 in [JSJ 
given in Section 3 and 4 therein, except that the proof of Lemma 4 of |JS] has to be refined 
in order to deal with strong mixing (jlj) instead of independent potential values V(S n uj). The 
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conclusion of the following lemma is hence exactly the same as of Lemma 4 of [JSj . and we 
thereby consider the proof of Theorem [2] to be complete. 

Let us set U = {E G C | E < $le(E) < E\ , \5sm(E)\ < 1}. Furthermore introduce the 
transfer matrices over several sites: 



k-l 



E 

\,S n UJ 5 



k > m 



Furthermore, T^ w (k,m) = (T^ w (m,k)) 1 for k < m and 7^(m, m) = 1. 
Lemma 5 Let E G U and N G N. Then there is a constant C such that the set 

Q N (E) = { uj e ft 

satisfies 



max IIT*>,1)|| 2 > e CN " 



0<n<N 



P(Q N {E)) > 1 - e 



-CN? 



Proof. For sake of notational symplicity, we will drop the index A on the transfer matrices 
Let us fix E G U and N G N and then split N into ^ pieces of length N 3 = N + N 1 + 2N 2 . 



For j — 0, . . . , j^, we consider the following events: 



u e Q 
uj g f2 



T*(jN 3 + N ,jN 3 )\\ < e^° N ° 
T*(jN 3 + N ,jN 3 )\\ < e*» N ° 



First we note that uniformly in uj and for some 71 > 

\\T*{n,m)\\ < e 7l '"~ m ' . 

Therefore the hypothesis (fl2|) implies as in the proof of Lemma 3 of [JSJ that, for E G U and 
iV G N, we have 



P(fi ) < 1-po < 1 



Po > 



(43) 



To shorten notations let us define ttj = 7ryN 3 -N2,jN 3 +N +N2] an d T^j = Xfi^j + N ,N 3 j). 
Using the quasi- locallity of g(u) = we get 



jN 3 +N -l / l-l 
l=jN 3 \k=jN 3 



q-E _ rr-E 

2 S l UJ 1 S l 1T J (uj) 



fjN 3 +N -l 

n r s k *M 

k=l+l 
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where C = C(g) as in Now choosing N2 = cNq for an adequate constant c, it follows that 



\T E -T E II < P h° N ° 



Therefore for a; € 



||r/ Hj || < ||7^|| + e^'° N ° < 2e5 7oAr ° < e ^ loN ° 

for ATq large enough, implying VP- C By a similar calculation, one obtains the second 
inclusion of 

n; c nj c nj. (44) 

By (1431) this implies 

P(fiJ) < P(J$ < 

Now clearly is 7^- = 7ry7v 3 _Ar 2j Ar 3+ jv 0+A r 2 )]-measurable. Therefore the strong mixing condition 
© implies that P(fiJ n Q\) < P(fiJ)P(fiJ) (1 + CiVf a ) < (1 - p ) 2 (1 + CWf a ). At the next 
step, one obtains P(fiJ n Q\ H < (1 - Po) 3 (1 + CN{ a ) 2 . Iteration and (ggj therefore give 

p( n n °) ^ p ( n ^) ^ ((i-po)(i+^r a ))^ . 

j=0,...,N/N 3 j=0,...,N/N 3 

Now let us choose N% sufficiently large such that 1 — p\ = (1 — po)(l + CN^ a ) < 1. Then 



G f2 



o< max^ ||7^(jW 3 + Ar ,jW 3 )|| 2 < e*" j. ) < 



Furthermore T E (jN 3 + N ,jN 3 ) — T®(jN s + N Q , l)T®(jN s , As A = EC implies either 
> or ||C|| — ll^ll 5 f° r arbitrary matrices, and = for A eSL(2,C), it 

therefore follows that 



uj e Q 



max max{||Tj ? (jiV3,l)|| 2 ,||Tj ? (jiV3 + iV ,l)|| 2 } > e 

0<j<N/N 3 



N i 

is greater or equal than 1 — (1 — pi) N z . Choosing N = cN? with adequate c concludes the 
proof. □ 
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